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Abstract 

It is commonly held that a necessary condition for the existence of solitons in 
nonlinear-wave systems is that the soliton's frequency (spatial or temporal) 
must not fall into the continuous spectrum of radiation modes. However, this is 
not always true. We present a new class of codimension-one solitons (i.e., those 
existing at isolated frequency values) that are embedded into the continuous 
spectrum. This is possible if the spectrum of the linearized system has (at 
least) two branches, one corresponding to exponentially localized solutions, 
and the other to radiation modes. An embedded soliton (ES) is obtained when 
the latter component exactly vanishes in the solitary-wave's tail. The paper 
contains both a survey of recent results obtained by the authors and some new 
results, the aim being to draw together several different mechanism underlying 
the existence of ESs. We also consider the distinctive property of semi- stability 
of ES, and moving ESs. Results are presented for four different physical models, 
including an extended 5th-order KdV equation describing surface waves in 
inviscid fluids, and three models from nonlinear optics. One of them pertains to 
a resonant Bragg grating in an optical fiber with a cubic nonlinearity, while two 
others describe second-harmonic generation (SHG) in the temporal or spatial 
domain (i.e., respectively, propagating pulses in nonlinear optical fibers, or 
stationary patterns in nonlinear planar waveguides). Special attention is paid 
to the SHG model in the temporal domain for of competing quadratic 

and cubic nonlinearities. An essential new result is that ES is, virtually, fully 
stable in the latter model in the case when both harmonics have anomalous 
dispersion. 
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1 Introduction 



Recent studies have revealed a novel type of solitary waves (which we will 
loosely called "solitons", without assuming integrability of the underlying 
models) that are embedded into the continuous spectrum, i.e., the soliton's 
internal frequency is in resonance with linear (radiation) waves. Generally, 
such a soliton should not exist, one finding instead a quasi- (delocalized) soli- 
ton with nonvanishing oscillatory tails (radiation component) [|TJ. Neverthe- 
less, bona fide (exponentially decaying) solitons can exist as codimension-one 
solutions if, at discrete values of the (quasi-) soliton's internal frequency, the 
amplitude of the tail exactly vanishes, while the soliton remains embedded into 
the continuous spectrum of the radiation modes. This requires the spectrum of 
the corresponding linearized system to consist of (at least) two branches, one 
corresponding to exponentially localized solutions, and the other to oscillatory 
radiation modes. In terms of the corresponding ordinary differential equations 
(ODEs) for the traveling-wave solutions, the origin must be a saddle- centre, 
that is its linearisation gives rise to both real and pure imaginary eigenvalues. 

Examples of such embedded solitons (ESs) were found in water-wave models 
p| and in several nonlinear-optical ones, including a Bragg-grating model with 
the wave-propagation (second-order-derivative) linear terms taken into regard 
||, and a model of the second-harmonic generation (SHG) in the presence of 
a Kerr nonlinearity The term "ES" was proposed in the latter work. 

ESs are interesting for several reasons, firstly because they frequently appear 
when higher-order (singular) perturbations are added to the system, which 
may completely change its soliton spectrum (see, e.g., 0). Secondly, optical 
ESs have considerable potential for applications, just because they are isolated 
solitons, rather than members of continuous families. Finally, and most cru- 
cial for physical applications, ESs are semi-stable objects. That is, as argued 
in Ref. [ffl analytically, in a general form applicable to ESs in any system, 
and checked numerically for the SHG model with the additional defocusing 
Kerr nonlinearity, ESs are stable in the linear approximation, but do have a 
slowly growing (sub-exponential) one-sided nonlinear instability (see Section 
| below). 

In the next section we discuss the existence and stability of ESs in four dif- 
ferent nonlinear partial-differential-equation (PDE) models. Mathematically 
speaking, in each case the reduced traveling-wave or steady-state ODEs have 
the structure of a fourth-order, reversible Hamiltonian system. In accord with 
what was said above, a parameter region we focus on is where the origin (trivial 
fixed point) in these ODE systems is a saddle-centre. That is, after diagonaliz- 
ing the system, one two-dimensional (2D) component of the dynamical system 
gives rise to imaginary eigenvalues ±iu (corresponding to a continuous radia- 
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tion branch in the linear spectrum of the PDE system), and the other to real 
eigenvalues ±A (corresponding to a gap in the radiation spectrum). 

All four models considered in this paper share the feature that ESs exist as 
non-generic, codimension-one solutions. Section 3 then presents three different 
views as to why this should be, which together provide for general insight 
into the existence and multiplicity of ESs. Section f| goes on to discuss their 
stability, arguing that ESs, in general, may be neutrally stable linearly, but 
suffer from a one-sided sub-exponential instability. This situation has been 
termed semi- stability in Q . Section [5] treats "moving" embedded solitons and, 
finally, Section (|6]) draws conclusions and briefly discusses potential physical 
applications of ESs in optical memory devices. 



2 Physical Examples 



2.1 An extended 5th-order KdV equation 



One of the first systems in which ESs were found (although without being given 
that name) is an extended 5th-order Korteweg - de Vries (KdV) equation 
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which with \x = reduces to the usual 5th-order KdV equation studied by a 
number of authors, see, e.g., Refs. |||7[]. The extended form ([I]) may be derived 
via a regular Hamiltonian perturbation theory from an exact Euler-equation 
formulation for water-waves with surface tension ||. Looking for traveling- 
wave solutions u(x — ct), integrating once, setting the constant of integration 
to be zero, and absorbing the linear term ~ u x by redefining a, one arrives at 
the following ODE (the prime stands for d/d(x — ct)), 
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(2) 



When a < 0, Eq. (|2]) is in the ES regime, since the linearization around the 
origin, u — 0, yields both real and imaginary eigenvalues. Note that, for the 
particular case \i — 0, it has been proved that there are no dynamical orbits 
homoclinic to the origin ||. Nevertheless, at least in the limit a — > —0, the 
system does possess a large family of homoclinic connections to periodic orbits, 
rather than to fixed points ||10|| . The latter generic solutions correspond to the 
above-mentioned delocalized solitary waves |TJ. 
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Fig. 1. The region of existence of the embedded solitons for Eq. (|2|) 

However, ESs do exist in the case \x = 1 ||, where one can find the explicit 
family of homoclinic-to-zero solutions, 



u (t) = 3 ( b + I) sech 2 ( \ 3{2 \ +1) t) , a= |(26 + 1)(6 - 2), 6 > -1/2.(3) 



2/ V 4 ' 5 



Moreover, this curve (family) of ESs was found in || to be only the first in a 
countable set of curves that appear to bifurcate from a = at a discrete set 
of negative values of b, see Fig. |l| Note that there are numerical difficulties in 
computing up to the limit point a = for b < 0; this is because, as will be 
motivated in Section |3]3] below, the bifurcation of these solutions from a = 
is a 'beyond-all-orders' effect 



We remark that Grimshaw and Cook IJTTJ found a similar bifurcation-type 
phenomenon in a system of two coupled KdV equations, although they did 
not explicitly identify the vanishing of the tail amplitude of the delocalized 
solitary waves as what we now call ESs. Also, Fujioka and Espinosa Jl2[ found 
a single explicit ES in a higher-order NLS equation with a quintic nonlinearity. 



Note also the following feature of this family of ES states. The first member 
of the family, viz., the explicit solution (|3|)) is a fundamental "ground state", 
while all other solution branches represent ' "excited states" , having the form 
of the ground-state soliton with small-amplitude "ripples" superimposed on 
it. As yet unpublished numerical results suggest that the ground-state soliton 
appears to be dynamically stable (in the sense described in Section f| below). 
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2.2 A generalized Massive Thirring model 



The first example of ESs in nonlinear optics was found in a generalized Thirring 
model (GTM), introduced long ago in 0,|T1]. ESs appear in this model when 



additional wave-propagation (second-derivative) terms are included (k being 
the carrier wavenumber), so that the model takes the form 



iu t + iu x + (2A;)- 1 {u xx - u tt ) + [(l/2)|u| 2 + M 2 ] u + v = 0, 
iv t - iv x + (2k)- 1 (v xx - v tt ) + [(1/2)|^| 2 + \u\ 2 ] v + u = 0. 



(4) 



Here u(x,t) and v(x,t) are right- and left- traveling waves coupled by reso- 
nant reflections on the grating. Soliton solutions are sought for as u(x,t) = 
exp(-iAut) {/(£) 5 v(x,t) = exp(— iAuit) V(0, where £ = x — ct, c and Auj 
being velocity and frequency shifts. The functions {7(0 and V(0 satisfy ODEs 

X U + z(l - C)U' + DU" + [(1/2)|{7| 2 + \V\ 2 } {7 + ^ = 0, 

( 5 ) 

X V - z(l + C)V + DV" + [(l/2)\V\ 2 + \U\ 2 ] V + U = 0, 



where x = Ac<j + (Au) 2 /2k, the effective velocity is C = (1 + Au/k)c, and an 
effective dispersion coefficient is = (1 — c 2 ) /2fc. The same equations were 
derived in for: (i) temporal-soliton propagation in nonlinear fiber gratings, 
including spatial-dispersion effects; and (ii) spatial solitons in a planar waveg- 
uide with a Bragg grating in the form of parallel scores, taking diffraction into 
regard, with t replaced by the propagation coordinate z, while x is the trans- 
verse coordinate. If we set C = then Eqs. (|5|) admit the further invariant 
reduction U = V*, leading to a single ODE ||, 

DU" + iU' + X U + (3/2)\U\ 2 U + U* = 0, (6) 



which is equivalent to a (Hamiltonian and reversible) system of four first-order 
equations for the real and imaginary parts of U and U' . In this system, the 
origin is a saddle-centre, provided that D > and \x\ < 1- In it was found 
numerically that Eq. (§) admits exactly three branches of the fundamental ESs 
(in contrast to Eq. (p where countably many ES states have been found). Also 
moving (c ^ 0) ESs, satisfying the 8th-order system (|5|), have also recently 
been found in Ref . [Hj] ; see Section [5] below for more details. 
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Fig. 2. (a) The (k,q) parameter plane of the three-wave model (0). The linear 
analysis (summarized in the inset boxes) shows that ESs can occur only in the 
region between the bold lines. The bundle of curves emanating from the point 
(k = 1, q = —4) are branches of embedded-soliton solutions with c = 0. The panels 
(b)-(e) depict solutions at the labeled points. 

2.3 A three-wave interaction model 



ESs can be found in far greater abundance in a model for spatial solitons, 
assuming a planar waveguide with a quadratic (x ) nonlinearity |l6fl , where 
two fundamental-harmonic (FH) waves vi t2 are coupled by the Bragg reflec- 
tions from a set of parallel scores. These two waves then interact nonlinearly, 
and generate a third wave, the second-harmonic (SH), with its wave- vector 
equal to the sum of those of the two FH components. The set of equations 
are: 

i(v lj2 )z ± i(vi, 2 )x + v 2 ,i + v 3 vl 1 = 0, 
2i(v 3 ) z - qv 3 + D(v 3 ) xx + v x v 2 = 0. 



Here v 3 is the SH field, x is the normalized transverse coordinate, q is a 
mismatch parameter, and D is an effective diffraction coefficient. Solutions 
to Eq. (0) are sought in the form Vi >2 (x, z) = exp(ikz)ui >2 (^), v 3 (x,z) = 
exp(2ikz) u 3 , with £ = x — cz, c being the slope of the soliton's axis relative 
to the propagation direction z. In Ref. ||17|| , many ESs of the zero-walkoff type 
(c = 0) were found (summarized in Fig. |2|), in which case the ODEs reduce, as 
in Eqs. (0), to a fourth-order real system. When c ^ 0, one cannot assume all 
the amplitudes to be real. In this case, one finds "moving" ESs as solutions 
to an eighth-order real ODE system (see Section below). 
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2.4 A second-harmonic-generation system 



The model we shall study in most detail is that for which the term 'ES' was 
first proposed in viz., a nonlinear optical medium with competing quadratic 
and cubic nonlinearities [HSIIS 



iu z + {l/2)u tt + u*v + 7i(|u| 2 + 2\v\ 2 )u = 0, 
iv z - (l/2)5v tt + qv + (l/2)u 2 + 2 72 (H 2 + 2\u\ 2 )v = 0. 



Here, u and v are FH and SH amplitudes, —5 is a relative SH/FH dispersion 
coefficient, q is a phase- velocity mismatch, and 7^2 are cubic (Kerr) nonlinear 
coefficients. In the absence of the Kerr nonlinearities, these equations are the 
same as those use by Karamzin and Sukhorukov in 1974 to obtain their 
famous soliton solution, in which both the FH and SH fields are propor- 
tional to sech 2 . A detailed analysis of the higher-order soliton solutions in that 



model with purely quadratic nonlinearity can be found in Ref. pi| . However 
the solutions that we will consider here are in a different class, in that the FH 
field will be more like a sech than a sech 2 . 

The particular case of Eqs. (H) with 5 = —1/2 is specially important, as 
it corresponds, with t replaced by the transverse coordinate x, to a second- 
harmonic-generation model in the spatial domain (in fact, in a nonlinear planar 
optical waveguide). In this special case, the model (^) is Galilean invariant, 
which allows one to generate a whole family of "moving" solitons from the 
single zero-walkoff one [22~]. At all other values of 5, construction of a "moving" 



(nonzero- walkoff) soliton is a nontrivial problem. 

Stationary solutions to Eq. (§) are sought for in the form u = U(t) exp(ikz), 
v = V(t) exp(2ikz), where k is real, and U, V satisfy ODEs 

(1/2)17" -W + U*V + n{\U\ 2 + 2\V\ 2 )U = 0, 
-(1/2)5V" + (q - 2k)V + (l/2)77 2 + 2 72 (|Vf + 2\U\ 2 )V = 0. 



In Ref. U, ES solutions to these equations were found for 5 > and 7 12 < 
(which implies anomalous and normal dispersions, respectively, at FH and 
SH, and self-defocusing Kerr nonlinearity. Alternatively, the same case may 
be physically realized as the normal and anomalous dispersions at FH and SH 
and self-focusing Kerr nonlinearity). In the same work, stability of these ESs 
was studied in detail. 

Here we shall present new results for 5 < 0, which corresponds to a more 
common case where the dispersion has the same sign at both harmonics. The 
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0.5 g 0.0 



Fig. 3. Two branches, in the (5, q) plane, of ES solutions to the second-harmonic 
generation model (||) with k = 0.3 and 71 = 72 = 0.05. The insets depict the profiles 
of U(t) (positive component) and V(t) (negative one). 

results can be naturally displayed in the form of curves of ESs in the (S, q) or 
(q, k) parameter planes, see Figs. [3] and [| below. 



3 Existence 



We now give three distinct explanations of why and how an ES may exist. 



3. 1 Nonlinearizablity 



Consider first the ODE system (|9]). It is important to note that the system 
gets fully decoupled in the linear approximation, the linearization of its sec- 
ond equation immediately telling one that no true soliton (with exponentially 
decaying tails) can exist inside the continuous (radiation- mode) SH spectrum. 
However, it may happen that the tail of the soliton's SH component decays 
at the same rate as the square of the tail of the FH component. In that case, 
the second equation of the system ([5|) is nonlinearizable, which opens the way 
for the existence of truly localized solitons inside the continuous spectrum. 
Note from the profiles of the solutions in Figs. ^ and ^, that the ^/-component 
appears to decay to zero much faster than U, in accordance with this nonlin- 
earizability property. 

Let us remark on some qualitative features of the ES branches presented in 
Fig. |^. First, the computed branches appear to end in "mid air". At the low q 
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Fig. 4. The branch of 6 = —0.5 ESs in the (£;,g)-plane corresponding to the 
right-hand branch in Fig. |3| (with 71 = 72 = 0.05), the insets similarly showing 
the soliton's profile. 

end, we have a boundary. The condition for the origin of the ODE system (H) 
to be a saddle-centre is k < q/2. For q below this limit, one may linearize (^), 
and find regular (non-embedded) solitons, for any q, provided k > 0. These 
latter solitons continuously match into the embedded solitons at q = k/2. 
Thus, since we took k = 0.3, we have that both branches of ESs will end at 
q = 0.6. At the high q end, there were numerical difficulties in continuing the 
branches, using the software AUTO [p^j , for precisely the same reasons why 
the computations presented in Fig. |l]were difficult at a = 0~ for b < 0. As one 
can see in the figure, one does have the amplitudes of both U and V tending 
to zero (note also that the V^-component is much smaller, in accordance with 
the nonlinearizability principle). 

Second, note that for high q, both ES solutions appear to be single-humped 
and fundamental. But upon going to smaller q, where the amplitudes of both 
U and V become larger, it is apparent that the right-hand branch is a higher- 
order state, with a superimposed ripple in the V component (akin to the 
second branch in Figs. |I| and |2|). 

Finally, we remark that only this second branch passes through the physically 
significant value of 5 = —1/2. For the spatial ESs, corresponding to 5 = —1/2, 
one is interested in how the solutions will vary as a function of k as well. That 
is shown in Fig. |J Again, the low-g limit corresponds to the boundary of the 
saddle-centre region is k = q/2. 
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3.2 Homoclinic orbits to saddle-centres 



Next, let us try to understand why ESs should be of codimension-one. We start 
from a general fourth-order Hamiltonian and reversible (invariant with respect 
to reflections of time and 'velocity' variables) ODE system. Linearizing about 
the zero solution (the origin), we assume the system to have a saddle-centre 
equilibrium. Thus, we have 



x = f(x), x e M\ 3R, R 2 = id, Rf(Rx) = -f(x), (10) 
3H(x) = const, along solutions, eigs(_D/(0)) = {A, — A, iu, — iu}. 

For example, for the ODE (0) the reversibility operator R is given by 

(it, u , u' , u") i— > (it, —u, u' , —u), 

and for the system @ by 

R:(U,U',V,V')^(U,-U',V,-V'). 

A homoclinic orbit connecting such a saddle-centre equilibrium to itself is 
formed by a trajectory simultaneously belonging to the one-dimensional un- 
stable and stable manifolds of the origin. Both of these manifolds lie in a 3D 
phase space H(0). Therefore, were it not for the reversibility, such homoclinic 
orbits would be of codimension-two in general |24j| , since we require the coin- 



cidence of two lines in the three-dimensional space. But, reversible homoclinic 
solutions (i.e., solutions that somewhere intersect the fixed-point set of the 
reversibility) are of codimension-one. This is because the unstable manifold 
and fix(i?) n H(0) are both one-dimensional, and we only require a point in- 
tersection between them. Hence, varying two parameters, we should expect to 
see ESs occurring along lines in the corresponding two-dimensional parameter 
plane. Moreover, the solutions themselves must be reversible; asymmetric ESs 
would be of a higher codimension still. 

Mielke, Holmes and O'Reilly proved a general theorem valid in the neigh- 
borhood of such a curve in the parameter plane of reversible saddle-centre ho- 
moclinic orbits: under a sign condition, essentially governing how reversibility 
and the Hamiltonian interact, they showed that there will be an accumulation 
of infinitely many curves of iV-pulse "bound states" of the primary homoclinic 
orbit, for each N > 1. Note that, for systems that are reversible and also have 
odd symmetry (such as @), the sign condition is always satisfied by virtue 
of the system's admitting both reversibilities R and —R (i.e., the model (|6]) 
is symmetric too under U — > —U, and hence is also invariant under the re- 
versibilities (U,U') -> (U*,-U'*) and (U,U') -> (-U*,U')). Here, the sign 
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condition determines whether there are "up-up" or "up-down" bound states. 
In Ref. [U, a large number of the bound states of the "up-down" type were 
found for the generalized massive Thirring model (|6|) in agreement with this 
theory. We also remark that Buryak and co-workers (see [p6fl ) found a simi- 
lar discrete sequence of bound states of 'nonexistent' dark solitons in a SHG 
model and a higher-order NLS equation, however none of these solutions were 
linearly stable. 

So far, in every case that we have investigated, we have found the higher- 
order ESs to be unstable against linear perturbations. A distinction should be 
stressed between what one would call "bound-states" — which are like several 
copies of a fundamental soliton placed end to end — and the higher-order 
solitons, such as those displayed in Figs. [I] and [| which are like a fundamental 
with internal ripples. At the moment, there seems to be no connection between 
these states, although it is still may happen that, as some parameter is varied, 
a continuous branch may connect solutions of the two different types. 



3.3 A singular limit 



We now look at a mechanism which explains how fundamental ESs (possibly 
with ripples) may appear from the singular limit A —>■ + . Such a limit for the 
general class of systems ([!(]) has been studied using the normal-form theory 
by Lombardi |27H 29[1 , incorporating careful estimation of various exponentially 
small terms (cf. related results obtained using exponential asymptotics, e.g. 



T0| , pTT|1 ) . A crucial additional ingredient we shall add to the Lombardi's work is 
that A and u are assumed to play the role of two independent parameters. We 
provide here only an oversimplified sketch, more details will appear elsewhere 
01. 



The appropriate normal form is [27] 



Xi = x 2 , 

x 2 = X 2 xi - {3/2)xj - b x {x\ + xl) +pN 2 {x; A), 
x 3 = —Xi{oj + b 2 xi) +pN 3 (x; A), 

x 4 = x 3 (uj + b 2 xi) +pN 4 (x; A) 



11^ 



where &i j2 are w-dependent constants to be determined for a particular sys- 
tem, and Ni are higher-order (remainder) terms that break (for nonzero p) 
the completely-integrable structure of the truncated normal form. It is not 
difficult to see that the truncated system possess a sech 2 -like homoclinic con- 
nection to the origin, whose amplitude is 0(A). The key question is whether 



11 



this homoclinic orbit persists under the inclusion of the remainder terms. This 
question can be posed in terms of the vanishing of a certain Melnikov inte- 
gral (whose vanishing measures the splitting distance between the stable and 



unstable manifolds), which, after a lengthy calculation f28|, can be written as 

P_ 

X- 



^ exp(-^7r/A) (A(AT 3) iV 4 , oj)(1 + 0(A) + 0(p)), 



where A(N 3 , N^,uj) can be computed explicitly for each monomial which is 
pure in X\ and x 2 in the Taylor-series expansion of either N 3 or N4. 

Now, something beautiful happens because, for each such monomial, A turns 



out to be [28| a (single- signed (^-dependent) constant multiple of either a 



Bessel function, A ~ J n (A\^ujb~2) for b 2 > 0, or modified Bessel function A 



7 n (4yu;|&2|) for b 2 < 0, for some integer n. Recall the basic properties of 
the Bessel functions, according to which J n (x) has infinitely many zeros for 
x > 0, whereas the modified Bessel function I n (x) = i~ n J n (ix) has no zeros. 
Hence, notice the crucial role played by the coefficient b 2 (in the truncated, 
scaled normal form ([TT]) ) in the case when the remainder N consists of a single 
monomial in {x\,x 2 ). If b 2 > 0, there will be an infinite number of ui values 
corresponding to zeros of A and hence homoclinic solutions will exist for small 
p and A. However, if 62 < 0, then A is strictly non-zero and hence there are 
no homoclinic solutions to the origin. 

The coefficient 62 is easy to calculate in the examples with the quadratic 
nonlinearities, such as the extended 5th-order KdV model (0) with p = 1 
(see |30| for details). There it is found b 2 > and this entirely explains the 



approximately periodic sequence of points on the negative 6-axis at which 
the homoclinic solutions bifurcate from a = with the zero amplitude, as 
shown in Fig. |]. In contrast, for p = — 1, b% is negative and no ES bifurcates 
from a = 0. For models with purely cubic nonlinearity, such as the generalized 
Thirring model (g), it can be shown that all quadratic coefficients in the normal 
form ( IT ) vanish, hence a new, odd-symmetric normal form would need to be 



studied. This is left for future work. We mention also that for the SHG model 
(|^), although there are quadratic terms in it, the coefficient b 2 is identically 
zero, so this analytic technique gives no information. 



4 Stability 



The stability of the embedded solitons in the SHG model has been studied 
both numerically and analytically in 1J. It was shown that the fundamental 
(single-hump) ES is linearly stable, but nonlinearly semi-unstable, while all 
the multi-humped ES are linearly unstable. In the semi- stability analysis of 
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Fig. 5. Evolution of a perturbed soliton with positive energy perturbation 



(a i = 0.1,a2 = —0.1 in (12)) for (JS|) from the left-hand fundamental branch of 
Fig. | 



the fundamental ES, a crucial role is played by the energy, as ESs are isolated 
solutions with uniquely determined values of the energy, with the adjacent de- 
localized soliton states, on either side, having an infinite energy. It was argued 
in 0j (see also a similar argument given by Buryak p6| ) that perturbations 
which slightly increases the fundamental ES's energy can be safe, while pertur- 
bation which decreases the energy inevitably triggers a slow (sub-exponential) 
decay of the soliton into radiation. Thus, the weak instability of an ES is 
one-sided. This fact also follows from a simple argument that the usual expo- 
nential instability is always dual-sided. Equivalently, the usual instability is 
linear, while the weak one-sided instability of ESs must be nonlinear. 



The situation is the same in the present case where 5 < 0. A study of the 
linearized equation around the ES shows that the fundamental ES branch 
(the left one in Fig. 3) is linearly stable, while the branch of multi-humped ES 
(the right one in Fig. 3) is linearly unstable. The semi-stability argument for 
the fundamental ES branch also applies here. Positive energy perturbations 
can be safe, while negative energy perturbations trigger decay of ES. However, 
as we shall see below, for 5 < 0, this decay seems to be significantly slower 
than that found in Ref. || for 5 > 0. As was done there, we numerically 
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simulated the system (§), with the initial data 

u(0, t) = U{t) + «isech2t, v(0, t) = V(t) + a 2 sech2t , (12) 



where U(t) and V(t) is an ES solution on the left-hand branch of Fig. |3] (the 
values are 5 = -2.9292, q = 6.0556, k = 0.2936251 and oti = a 2 = 0.05). Figs. 
|5] and |G] depict, respectively, the effects of the positive-energy and negative- 
energy perturbations. In both cases, we observe fast oscillations on top of 
slow ones in the evolution of \u\ and \v\ at t = 0. These oscillations are very 



similar to those reported in []22| for perturbed solitons in the standard SHG 
system with no cubic nonlinearity (note that those solitons are ordinary ones, 
rather than ES). In that case, the fast oscillations were attributed to an intrin- 
sic mode, while the slow oscillations were attributed to beatings between the 
intrinsic mode and a quasi-mode. The latter one is localized in the FH compo- 
nent, but resonates with the continuous spectrum in the SH component. Both 
oscillations could last for a very long time, even though they were expected 
to eventually decay due to a very weak radiation damping. 

We believe that similar mechanisms are also at work in our model. However, 
there are important differences because of the fact that the solitons in our 
model are embedded, and those in the model considered in Ref. |22| were not. 



On the other hand, there are important similarities because perturbations of 
an ES naturally resonate with the continuous spectrum of the SH compo- 
nent, and we do see such slow oscillations as well. Although according to the 
semi-stability argument, a negative perturbation of an ES, as shown in Fig. ^ 
would eventually decay, it is a very slow decay. Detailed examination of the 
numerical solutions shows that the central pulse (the v component) in Fig. ^ 
keeps shedding oscillating tails into the far field. However, the tail amplitudes 
are extremely small (about 0.001 or smaller). This is why the expected decay 
of the perturbed ES is not obvious in that figure. Because of this, the actually 
observed evolution is dominated by the beating and internal oscillations, just 
as in Fig. |5| with a positive perturbation, and as in Ref. [22]. It will take an 



extremely long time for the pulse in Fig. [] to show considerable decay. In 
fact, it is clear that the ESs in the present model with 5 < are virtually 
stable, when compared to the previously considered case 5 > 0, where the 
semi-instability was a really observed feature. The relative stability of an ES 
for 5 < is a new result reported in the present paper, and its importance is 
quite obvious. 

If the ES were linearly unstable, the semi-stability argument would not apply. 
This is the case for the solutions investigated in |[26|| , and also for the right-hand 
branch depicted in Fig. || which corresponds to the case of a nonlinear planar 
optical waveguide in the spatial domain. To verify this, we have performed a 
time integration of the PDE, the results of which are presented in Fig. [7|. One 
can see the onset of a violent exponential instability. We have verified that 
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Evolution of a perturbed soliton with a negative energy perturbation 
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Fig. 7. The linear instability of a perturbed soliton belonging to the right-hand 
branch in Fig. g at S = -0.5, 7l = 72 = 0.05, q = 10 and k = 1.9298 

this solution does have exponentially unstable eigenmodes, from a numerical 
study of the linearization of Eq. (||) about this solution. 

In view of these results we conjecture the following. Fundamental ESs are, 
in general, linearly neutrally stable but semi-stable nonlinearly. The higher- 
order ESs (which usually have internal ripples in their profiles) are, generally, 
linearly unstable. Preliminary results for the extended 5th-order KdV ([I]) in- 
dicate qualitatively the same properties. This would also accord with previous 
numerical results for higher-order NLS equations that an isolated fundamental 
ES |12| is semi-stable whereas multihumped bound-states J26| are exponen- 
tially unstable. 
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Fig. 8. A branch of moving ESs bifurcating from one of the three simplest branches of 
the fundamental quiescent ESs in Eqs. (|^) with a = 1/2: the velocity (a), momentum 
(b), and energy (c) vs. D, with insets showing the shape of the moving solitons at 
labeled points. 

To conclude this section, we notice that multihumped solitons (which are not a 
subject of the present paper) of the ordinary (nonembedded) type were found 
in many models, see, e.g., early works |pl| , |32| dealing with the Langmuir waves 
in plasmas, and a recent work on the resonant three-wave interactions 



5 Moving embedded solitons 

We will now discuss a possibility that ESs in the generalized Massive Thirring 
model (HD may be moving at a non-zero velocity c (this discussion applies also 
to ESs in the three- wave system (0) with a non-zero walkoff). In this case, 
the reduction of the 8th-order ODE to a 4th-order one is no longer possi- 
ble. Moreover, since the 8th-order system is obtained by separating the real 
and imaginary parts of a 4th-order complex system, the spectrum will have 
a double degeneracy when c = 0. This means that, in the parameter region 
of interest in the (x, c, _D)-space, the linearization yields four pure imaginary 
eigenvalues, plus two with positive real parts and two with negative ones. 
A similar counting argument, as in Section ( |3.2|) , shows that reversible ho- 
moclinic orbits to such equilibria are of codimension two. Hence ESs lie on 
one-dimensional curves in the three-parameter space. Alternatively, this prop- 
erty can be explained as follows: in addition to the energy E, the full system 
also preserves the momentum P, so we can view a moving ES as being isolated 
in both invariants, i.e., the ES solution family is described by curves E(D) and 
P{D). Finally, we find that such curves can be found naturally as bifurcation 
points at c = from curves of the zero-velocity ESs. 
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Fig. 9. Two branches of "walking" (c 7^ 0) embedded solitons bifurcating from the 
ground-state c = branch (solid horizonal line) for the three- wave model ([?]): (a) 
the walkoff c, and (b) the energy flux E vs. the mismatch q. The inset in (b) shows 
that the two branches meet on the left and merge in a typical tangent bifurcation. 

Fig. ^| shows one such bifurcation occurring from one of the three branches 
of fundamental ESs for the Thirring model. The ES branches terminate (as 
at point "4") and beyond that, ES become ordinary (non-embedded) solitons 
(i.e., where the saddle-centre equilibrium becomes a pure saddle with only real 
eigenvalues) . 

Similarly, Fig. |S] shows two branches of moving ESs bifurcating from the 
ground-state ES (the nearly vertical line) in the three-wave model. Note that 
the two bifurcating curves become globally connected on the left, through a 
regular tangent (fold) bifurcation. 

Finally, let us turn to the second-harmonic model (§). As it was mentioned 
above, precisely at the value 5 = —1/2, this model gives rise to moving solitons 
in a trivial way, via the Galilean transformation. Search for moving ES in this 
model with 5 7^ —1/2 is the subject of ongoing work. 



6 Conclusion 



In this paper, we have given a brief overview of recent results that establish 
the existence of isolated (codimension-one) solitons embedded into the con- 
tinuous spectrum of radiation modes. A necessary condition for the existence 
of the embedded solitons is the presence of (at least) two different branches in 
the spectrum of the corresponding linearized system, so that one branch can 
correspond to purely imaginary eigenvalues, and another to purely real ones. 
The fundamental (single-humped) embedded solitons are always stable in the 
linear approximation, being subject to a weak sub-exponential one-sided insta- 
bility. Moving embedded solitons may also exist as codimension-two solutions. 
Moreover, bound states in the form of multi-humped embedded solitons exist 
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too, but they are linearly unstable. 



We should remark that the work presented and reviewed here does not nec- 
essarily represent the very first time that the existence of ESs have been es- 
tablished in any physical model. For example, as stated earlier, multi-humped 
dark ESs were already observed in a SHG model and for higher-order NLS 
equations [^6| where they were found to be unstable. Fundamental ESs were 
also found for the latter with competing cubic and quintic nonlinearities 



and their existence was also suggested for a coupled KdV system |TT[, without 



evidence to suggest their stability. Also recent work |33|,|34| has shown the exis- 
tence of front solutions to higher-order Frekel-Kontorova models which are in 
resonance with the linear spectrum; so called 'embedded kinks'. Nevertheless, 
we believe the results presented in this paper — the new stability results for 
fundamental ESs in Section 4, together with the three distinct views given in 
Section 3 as to why ESs should exist, and the similarities found in Section 2 
between their existence properties for four distinct models — provide a new 
theory for ESs phenomenon in their own right. 



Moreover, the very fact that ESs are isolated states suggest that they may 
find potential application in photonics, such as in all-optical switching. For 
example, taking the examples in Sections 2.2 and 2.3, switching from one ES 
state to a neighboring one with a smaller energy might be easily initiated by 
a small perturbation, in view of the semi-stability inherent to ESs. Switching 
between two branches of moving ESs with c 7^ might be quite easy to realize 
too, due to the small energy-flux and walkoff differences between them. There 
remains much work to be done in investigating these potential aplications 
further. 

Many issues concerning the embedded solitons remain open and are a subject 
of ongoing investigations. In particular, immediate questions arise concerning 
moving ESs, and interactions caused by collisions between them. 
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supported in part by the NSF and the AFOSR. The research of DJK was 
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through the University of Bristol. ARC is suppoted by an Advanced Fellowship 
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